Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lysozyme”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Fluid Physics and Macromolecular Crystal Growth in Microgravity

The molecular structure of biological macromolecules is important in understanding how these molecules work and has direct application to rational drug design for new medicines and for the improvement and development of industrial enzymes. In order to obtain the molecular structure, large, well formed, single macromolecule crystals are required. The growth of macromolecule crystals is a difficult task and is often hampered on the ground by fluid flows that result from the interaction of gravity with the crystal growth process. One such effect is the bulk movement of the crystal through the fluid due to sedimentation. A second is buoyancy driven convection close to the crystal surface. On the ground the crystallization process itself induces both of these flows. Buoyancy driven convection results from density differences between the bulk solution and fluid close to the crystal surface which has been depleted of macromolecules due to crystal growth. Schlieren photograph of a growing lysozyme crystal illustrating a 'growth plume' resulting from buoyancy driven convection. Both sedimentation and buoyancy driven convection have a negative effect on crystal growth and microgravity is seen as a way to both greatly reduce sedimentation and provide greater stability for 'depletion zones' around growing crystals. Some current crystal growth hardware however such as those based on a vapor diffusion techniques, may also be introducing unwanted Marangoni convection which becomes more pronounced in microgravity. Negative effects of g-jitter on crystal growth have also been observed. To study the magnitude of fluid flows around growing crystals we have attached a number of different fluorescent probes to lysozyme molecules. At low concentrations, less than 40% of the total protein, the probes do not appear to effect the crystal growth process. By using these probes we expect to determine not only the effect of induced flows due to crystal growth hardware design but also hope to optimize crystallization hardware so that destructive flows are minimized both on the ground and in microgravity.

Pusey, M.↗
NASA NTRS2001

Investigating the Effect of Impurities on Macromolecule Crystal Growth in Microgravity

Chicken egg-white lysozyme (CEWL) crystals were grown in microgravity and on the ground in the presence of various amounts of a naturally occurring lysozyme dimer impurity. No significant favorable differences in impurity incorporation between microgravity and ground crystal samples were observed. At low impurity concentration the microgravity crystals preferentially incorporated the dimer. The presence of the dimer in the crystallization solutions in microgravity reduced crystal size, increased mosaicity and reduced the signal to noise ratio of the X-ray data. Microgravity samples proved more sensitive to impurity. Accurate indexing of the reflections proved critical to the X-ray analysis. The largest crystals with the best X-ray diffraction properties were grown from pure solution in microgravity.

Snell, Edward H.↗
NASA NTRS2000

A Compact X-Ray System for Macromolecular Crystallography

We describe the design and performance of a high flux x-ray system for a macromolecular crystallography that combines a microfocus x-ray generator (40 micrometer full width at half maximum spot size at a power level of 46.5 W) and a collimating polycapillary optic. The Cu Ka lpha x-ray flux produced by this optimized system through a 500,um diam orifice is 7.0 times greater than the x-ray flux previously reported by Gubarev et al. [M. Gubarev et al., J. Appl. Crystallogr. 33, 882 (2000)]. The x-ray flux from the microfocus system is also 2.6 times higher than that produced by a rotating anode generator equipped with a graded multilayer monochromator (green optic, Osmic Inc. CMF24-48-Cu6) and 40% less than that produced by a rotating anode generator with the newest design of graded multilayer monochromator (blue optic, Osmic, Inc. CMF12-38-Cu6). Both rotating anode generators operate at a power level of 5000 W, dissipating more than 100 times the power of our microfocus x-ray system. Diffraction data collected from small test crystals are of high quality. For example, 42 540 reflections collected at ambient temperature from a lysozyme crystal yielded R(sub sym)=5.0% for data extending to 1.70 A, and 4.8% for the complete set of data to 1.85 A. The amplitudes of the observed reflections were used to calculate difference electron density maps that revealed positions of structurally important ions and water molecules in the crystal of lysozyme using the phases calculated from the protein model.

Gubarev, Mikhail↗
NASA NTRS2001

A Compact X-Ray System for Support of High Throughput Crystallography

Standard x-ray systems for crystallography rely on massive generators coupled with optics that guide X-ray beams onto the crystal sample. Optics for single-crystal diffractometry include total reflection mirrors, polycapillary optics or graded multilayer monochromators. The benefit of using polycapillary optic is that it can collect x-rays over tile greatest solid angle, and thus most efficiently, utilize the greatest portion of X-rays emitted from the Source, The x-ray generator has to have a small anode spot, and thus its size and power requirements can be substantially reduced We present the design and results from the first high flux x-ray system for crystallography that combine's a microfocus X-ray generator (40microns FWHM Spot size at a power of 45 W) and a collimating, polycapillary optic. Diffraction data collected from small test crystals with cell dimensions up to 160A (lysozyme and thaumatin) are of high quality. For example, diffraction data collected from a lysozyme crystal at RT yielded R=5.0% for data extending to 1.70A. We compare these results with measurements taken from standard crystallographic systems. Our current microfocus X-ray diffraction system is attractive for supporting crystal growth research in the standard crystallography laboratory as well as in remote, automated crystal growth laboratory. Its small volume, light-weight, and low power requirements are sufficient to have it installed in unique environments, i.e.. on-board International Space Station.

Ciszak, Ewa↗
NASA NTRS2002

Seeing the Heat: Preliminary Studies of Cryocrystallography Using Infrared Imaging

As preparation for an extensive study imaging the cryocooling process of macromolecular crystals we have demonstrated the ability to thermally image solid objects and liquids at temperatures far below 273 K. In the case of a large lysozyme crystal qualitative measurements show the cooling process to take about 0.6s with the cooling taking place in a wave from the face of the crystal nearest to the origin of the cryostream, to the point furthest away from the origin. Annealing of this lysozyme crystal, cooled under good cryoprotectant conditions, showed cold striations formed perpendicular to the cooling stream. These striations became more pronounced after successive annealing. Cryocooling of a non-cryoprotected crystal of glucose isomerase displayed an S-shaped cold front wave traveling across the sample. These preliminary results are qualitative but show the power of infrared imaging as a new tool for fundamental and practical cryocrystallography studies.

Snell, E. H.↗
NASA NTRS2002

A Proposed Model for Protein Crystal Nucleation and Growth

How does one take a molecule, strongly asymmetric in both shape and charge distribution, and assemble it into a crystal? We propose a model for the nucleation and crystal growth process for tetragonal lysozyme, based upon fluorescence, light, neutron, and X-ray scattering data, size exclusion chromatography experiments, dialysis kinetics, AFM, and modeling of growth rate data, from this and other laboratories. The first species formed is postulated to be a 'head to side' dimer. Through repeating associations involving the same intermolecular interactions this grows to a 4(sub 3) helix structure, that in turn serves as the basic unit for nucleation and subsequent crystal growth. High salt attenuates surface charges while promoting hydrophobic interactions. Symmetry facilitates subsequent helix-helix self-association. Assembly stability is enhanced when a four helix structure is obtained, with each bound to two neighbors. Only two unique interactions are required. The first are those for helix formation, where the dominant interaction is the intermolecular bridging anion. The second is the anti-parallel side-by-side helix-helix interaction, guided by alternating pairs of symmetry related salt bridges along each side. At this stage all eight unique positions of the P4(sub3)2(sub 1),2(sub 1) unit cell are filled. The process is one of a) attenuating the most strongly interacting groups, such that b) the molecules begin to self-associate in defined patterns, so that c) symmetry is obtained, which d) propagates as a growing crystal. Simple and conceptually obvious in hindsight, this tells much about what we are empirically doing when we crystallize macromolecules. By adjusting the growth parameters we are empirically balancing the intermolecular interactions, preferentially attenuating the dominant strong (for lysozyme the charged groups) while strengthening the lesser strong (hydrophobic) interactions. In the general case for proteins the lack of a singularly defined association pathway may lead to formation of multiple species, i.e., amorphous precipitation. Weak interactions, such as hydrogen bonds, are promiscuous, serving to strengthen rather than define specific interactions. Participation in an interaction sequesters that surface from subsequent interactions, and we expect the strongest bonds to form first. This model, its basis, how it fits into the currently understood osmotic second virial coefficient approach to crystallization, and what it suggests will be discussed.

Pusey, Marc↗
NASA NTRS2003

Perfection of Apoferritin Crystals: An Advanced X-Ray Imaging and Diffraction Study

Ferritin is a well-known iron-storage protein, and is a spherical shell that consists of 24 identical subunits packed in a 432 symmetry. The typically large protein size and its distinction from lysozyme as to chemical and physical characteristics make ferritin an attractive model protein for crystal growth and perfection investigation-as an alternative to the most widely studied lysozyme. In this contribution, the latest results obtained from coherence-based x-ray diffraction imaging and diffraction experiments will be presented on octahedral apoferritin (a demetalized form of ferritin) crystals grown from various growth conditions. Crystal specimens, which have the measured rocking-curve widths varying from a few arcseconds to several tens arcseconds (or more), are comparatively examined by intrinsically highly sensitive mapping of lattice perfection and defects. The richness of the observed defects and growth features offers insight into perfection and growth of protein crystals. Beautiful interference fringe patterns formed in diffraction images and fine oscillation structure of rocking curves observed will be discussed for understanding of physical origins and the underlying impact.

Hu, Z. W.↗
NASA NTRS2003

An Alternative Hypothesis for How Microgravity Improves Macromolecular Crystal Quality

There is a substantial body of experimental evidence, from this and other laboratories, that strongly suggests that for many proteins crystal nucleation and growth is by addition of associated species that are preformed by reversible concentration-driven self association processes in the bulk solution. We have developed a self-association model for the crystal nucleation and growth of the protein chicken egg lysozyme. The model accounts for the obtained crystal symmetry, explains the observed surface structures, and shows the importance of the symmetry obtained by self-association in solution to the process as a whole. This model also offers a possible mechanism for fluid flow effects on the growth process and how microgravity may affect it. While a single lysozyme molecule is relatively small an octamer in the 43 helix configuration (the proposed average sized growth unit) would have a M.W. approx. 115,000 and dimensions of 5.6 x 5.6 x 7.6 nm. Direct AFM measurements of growth unit incorporation indicate that units as wide as 11.2 nm and as long as 11.4 nm (a 24-mer) commonly attach to the crystal. AFM results from Weichmann et al. (Ultramicroscopy 86, 159-166, 2001) suggest that associated species of up to 40-mers in size add to the (101) faces. These measurements reflect the sizes of units that both added and desorbed from the crystal surface. The larger and less isotropic the associated species the more likely that it will be oriented to some degree in a flowing boundary layer, even at the low flow velocities measured about macromolecule crystals. On Earth, concentration gradient driven flow will maintain a high interfacial concentration, i.e., a high level (essentially that of the bulk solution) of solute association at the interface and higher growth rate. Higher growth rates mean an increased probability that misaligned growth units are trapped by subsequent growth layers before they can be desorbed and try again, or that the desorbing species is more likely to be smaller than the adsorbing species. In microgravity the extended diffusive boundary layer will lower the interfacial concentration. This results in a net dissociation of aggregated species that diffuse in from the bulk solution, i.e., smaller associated species, which are more likely able to make multiple attempts to correctly bind, yielding higher quality crystals.

Pusey, Marc↗
NASA NTRS2004

Solutal Convection Around Growing Protein Crystal and Diffusional Purification in Space

At least some protein crystals were found to preferentially trap microheterogeneous impurities. The latter are, for example, dimmer molecules of the crystallizing proteines (e.g. ferritin, lysozyme), or the regular molecules on which surfaces small molecules or ions are adsorbed (e.g. acetilated lysozyme) and modi@ molecular charge. Impurities may induce lattice defects and deteriorate structural resolution. Distribution of impurities between mother solution and gorwing crystal is defined by two interrelated distribution coefficients: kappa = rho(sup c, sub 2) and K = (rho(sup c, sub 2)/rho(sup c, sub 1)/rho(sub 2)/rho(sub 1). Here, rho(sub 2), rho(sub 1) and rho(sup c, sub 2) are densities of impurity (2) and regular protein (1) in solution at the growing interface and within the crystal ("c"). For the microheterogeneous impurities studied, K approx. = 2 - 4, so that kappa approx. - 10(exp 2) - 10(exp 3), since K = kappa (rho(sub 1)/rho(sup c, sub 1) and protein solubility ratio rho(sub 1)/rho(sub=p c, sub 2) much less than 1. Therefore, a crystal growing in absence of convection purifies mother solution around itself, grows cleaner and, probably, more perfect. If convection is present, the solution flow permanently brings new impurities to the crystal. This work theoretically addressed two subjects: 1) onset of convection, 2) distribution of impurities.

Lee, Chun P.↗
NASA NTRS2005

Liquid Between Macromolecules in Protein Crystals: Static Versus Dynamics

Protein crystals are so fragile that they often can not be handled by tweezers. Indeed, measurements of the Young modulus, E, of lysozyme crystals resulted in E approx. equals 0.1 - 1 GPa, the lower figures, 0.1 - 0.5 GPa, being obtained from triple point bending of as-grown and not cross-linked crystals sitting in solution. The bending strength was found to be approx.10(exp -2) E. On the other hand, ultrasound speed and Mandelstam-Raman-Brilloin light scattering experiments led to much higher figures, E approx. equals 2.7 GPa. The lower figures for E were found from static or low frequency crystal deformations measurements, while the higher moduli are based on high frequency lattice vibrations, 10(exp 7) - 10(exp 10) 1/s. The physical reason for the about an order of magnitude discrepancy is in different behavior of water filling space between protein molecules. At slow lattice deformation, the not-bound intermolecular water has enough time to flow from the compressed to expanded regions of the deformed crystal. At high deformation frequencies in the ultra- and hypersound waves, the water is confined in the intermolecular space and, on that scale, behaves like a solid, thus contributing to the elastic crystal moduli. In this case, the reciprocal crystal modulus is expected to be an average of the water protein and water compressibilities (reciprocal compressibilities): the bulk modulus for lysozyme is 26 GPa, for water it is 7 GPa. Anisotropy of the crystal moduli comes from intermolecular contacts within the lattice while the high frequency hardness comes from the bulk of protein molecules and water bulk moduli. These conclusions are based on the analysis of liquid flow in porous medium to be presented.

Chernov, A. A.↗
NASA NTRS2005

Growth from Solutions: Kink dynamics, Stoichiometry, Face Kinetics and stability in turbulent flow

1. Kink dynamics. The first segment of a polygomized dislocation spiral step measured by AFM demonstrates up to 60% scattering in the critical length l*- the length when the segment starts to propagate. On orthorhombic lysozyme, this length is shorter than that the observed interkink distance. Step energy from the critical segment length based on the Gibbs-Thomson law (GTL), l* = 20(omega)alpha/(Delta)mu is several times larger than the energy from 2D nucleation rate. Here o is tine building block specific voiume, a is the step riser specific free energy, Delta(mu) is the crystallization driving force. These new data support our earlier assumption that the classical Frenkel, Burton -Cabrera-Frank concept of the abundant kink supply by fluctuations is not applicable for strongly polygonized steps. Step rate measurements on brushite confirms that statement. This is the1D nucleation of kinks that control step propagation. The GTL is valid only if l* <Dk/vk, the diffusion path of a kink that has diffusivity Dk and average growth velocity vk. This is equivalent to supersaturations sigma less than approx. alpha/2l*, where alpha is the building block size. For lysozyme, sigma much less than (1%). Conventionally used interstep distance generated by screw dislocation, 19(omega)alpha/Delta(mu) should be replaced by the very different real one, approx.4l*. 2. Stoichiometry. Kink, and thus step and face rates of a non-Kossel complex molecular monocomponent or any binary, AB, lattice was found theoretically to be proportional to 1/(zeta(sup 1/2) + zeta(sup - 1/2)), where zeta = [B]/[A] is the stoichiometry ratio in solution. The velocities reach maxima at zeta = 1. AFM studies of step rates on CaOxalate monohydrate (kidney stones) from aqueous solution was found to obey the law mentioned above. Generalization for more complex lattice will be discussed. 3. Turbulence. In agreement with theory, high precision in-situ laser interferometry of the (101) KDP crystal face shows step bunching if solution flows parallel to the step flow. The bunch height increases with the distance the bunch travels, i.e. with the face size. However, when the flow rate, u, increases, at u greater than approx. 1 m / s , the average step bunch height decreases as 1/u. The pheonomenon is attributed to the turbulent rather than laminar viscous boundary layer where diffusivity Dt = 0.5u(sub tau),y, i.e. increases linearly with the distance y from the solid face. Friction velocity, u(sub tau) approx. u(sup 7/8). Dramatically larger rate of the mass/heat transport within the turbulent, as compared to the laminar, viscous layer will be discussed.

Chernov, A. A.↗
DOE OSTI2021

The exo-β- N -acetylmuramidase NamZ from Bacillus subtilis is the founding member of a family of exo-lytic peptidoglycan hexosaminidases

Endo-β-N-acetylmuramidases, commonly known as lysozymes, are well-characterized antimicrobial enzymes that catalyze an endo-lytic cleavage of peptidoglycan; i.e., they hydrolyze the β-1,4-glycosidic bonds connecting N-acetylmuramic acid (MurNAc) and N-acetylglucosamine (GlcNAc). In contrast, little is known about exo-β-N-acetylmuramidases, which catalyze an exo-lytic cleavage of β-1,4-MurNAc entities from the non-reducing ends of peptidoglycan chains. Such an enzyme was identified earlier in the bacterium Bacillus subtilis, but the corresponding gene has remained unknown so far. We now report that ybbC of B. subtilis, renamed namZ, encodes the reported exo-β-N-acetylmuramidase. A ΔnamZ mutant accumulated specific cell wall fragments and showed growth defects under starvation conditions, indicating a role of NamZ in cell wall turnover and recycling. Recombinant NamZ protein specifically hydrolyzed the artificial substrate para-nitrophenyl β-MurNAc and the peptidoglycan-derived disaccharide MurNAc-β-1,4-GlcNAc. Together with the exo-β-N-acetylglucosaminidase NagZ and the exo-muramoyl-L-alanine amidase AmiE, NamZ degraded intact peptidoglycan by sequential hydrolysis from the non-reducing ends. A structure model of NamZ, built on the basis of two crystal structures of putative orthologs from Bacteroides fragilis, revealed a two-domain structure including a Rossmann-fold-like domain that constitutes a unique glycosidase fold. Thus, NamZ, a member of the DUF1343 protein family of unknown function, is now classified as the founding member of a new family of glycosidases (CAZy GH171). NamZ-like peptidoglycan hexosaminidases are mainly present in the phylum Bacteroidetes and less frequently found in individual genomes within Firmicutes (Bacilli, Clostridia), Actinobacteria, and γ-proteobacteria.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Cellular Phosphate Sensing and Anion Binding by an Azacrown‐Calixpyrrole Hybrid

A hybrid receptor-sensor for anions originating from the merging of positively charged ammonium moieties for electrostatic attraction/stronger binding of azacrowns with directionality of calixpyrrole hydrogen bond donors for selectivity is investigated. As demonstrated this hybrid receptor-sensor shows a remarkable selectivity for orthophosphate even in the presence of other phosphates and anions found in cellular materials (K assoc H 2 PO 4 − >H 2 P 2 O 7 2− >AMP − ≫ADP 2− or ATP 3− over halides, nitrate, or hydrogen sulfate; all Na + salts in water) but also cellular polyphosphate or phospholipids. This selectivity is harnessed in a real-time monitoring of cell lysis by lysozyme, which releases orthophosphate and other phosphates and anions from the cells. This sensitive (LOD 0.4 μM) fluorescence-based microscale method compares favorably with the state-of-the-art techniques but can easily be practiced in a high-throughput screening (HTS) manner. The anion binding and selectivity in aqueous solutions were investigated by NMR and put in context with phosphate binding of the parent calix[4]pyrrole. The microscopic understanding of anion binding by the hybrid receptor was then obtained from a combination of density functional theory (DFT), classical molecular dynamics (MD) with explicit water solvation, and ab initio MD (AIMD) simulations. Correlating the NMR and fluorescence binding data with studies of solvation of the receptor, phosphate anion, and the resulting complex confirms the binding is largely driven by entropic component (TΔS) associated with receptor and anion desolvation.

Anions↗
DOE OSTI2023

A Highly Ordered Nitroxide Side Chain for Distance Mapping and Monitoring Slow Structural Fluctuations in Proteins

Abstract Site-directed spin labeling electron paramagnetic resonance (SDSL-EPR) is an established tool for exploring protein structure and dynamics. Although nitroxide side chains attached to a single cysteine via a disulfide linkage are commonly employed in SDSL-EPR, their internal flexibility complicates applications to monitor slow internal motions in proteins and to structure determination by distance mapping. Moreover, the labile disulfide linkage prohibits the use of reducing agents often needed for protein stability. To enable the application of SDSL-EPR to the measurement of slow internal dynamics, new spin labels with hindered internal motion are desired. Here, we introduce a highly ordered nitroxide side chain, designated R9, attached at a single cysteine residue via a non-reducible thioether linkage. The reaction to introduce R9 is highly selective for solvent-exposed cysteine residues. Structures of R9 at two helical sites in T4 Lysozyme were determined by X-ray crystallography and the mobility in helical sequences was characterized by EPR spectral lineshape analysis, Saturation Transfer EPR, and Saturation Recovery EPR. In addition, interspin distance measurements between pairs of R9 residues are reported. Collectively, all data indicate that R9 will be useful for monitoring slow internal structural fluctuations, and applications to distance mapping via dipolar spectroscopy and relaxation enhancement methods are anticipated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Intact mass spectrometry screening to optimize hydroxyl radical dose for protein footprinting

Hydroxyl radical protein footprinting (HRPF) using synchrotron radiation is a well-validated method to assess protein structure in the native solution state. In this method, X-ray radiolysis of water generates hydroxyl radicals that can react with solvent accessible side chains of proteins, with mass spectrometry used to detect the resulting labeled products. An ideal footprinting dose provides sufficient labeling to measure the structure but not so much as to influence the results. The optimization of hydroxyl radical dose is typically performed using an indirect Alexa488 fluorescence assay sensitive to hydroxyl radical concentration, but full evaluation of the experiment's outcome relies upon bottom-up liquid chromatography mass spectrometry (LC-MS) measurements to directly determine sites and extent of oxidative labeling at the peptide and protein level. A direct evaluation of the extent of labeling to provide direct and absolute measurements of dose and “safe” dose ranges in terms of, for example, average numbers of labels per protein, would provide immediate feedback on experimental outcomes prior to embarking on detailed LC-MS analyses. To this end, we describe an approach to integrate intact MS screening of labeled samples immediately following exposure, along with metrics to quantify the extent of observed labeling from the intact mass spectra. Intact MS results on the model protein lysozyme were evaluated in the context of Alexa488 assay results and a bottom-up LC-MS analysis of the same samples. In conclusion, this approach provides a basis for placing delivered hydroxyl radical dose metrics on firmer technical grounds for synchrotron X-ray footprinting of proteins, with explicit parameters to increase the likelihood of a productive experimental outcome. Further, the method directs approaches to provide absolute and direct dosimetry for all types of labeling for protein footprinting.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Persistent Protein Motions in a Rugged Energy Landscape Revealed by Normal Mode Ensemble Analysis

Proteins are allosteric machines that couple motions at distinct, often distant, sites to control biological function. Low-frequency structural vibrations are a mechanism of this long-distance connection and are often used computationally to predict correlations, but experimentally identifying the vibrations associated with specific motions has proved challenging. Spectroscopy is an ideal tool to explore these excitations, but measurements have been largely unable to identify important frequency bands. The result is at odds with some previous calculations and raises the question what methods could successfully characterize protein structural vibrations. Here we show the lack of spectral structure arises in part from the variations in protein structure as the protein samples the energy landscape. However, by averaging over the energy landscape as sampled using an aggregate 18.5 μs of all-atom molecular dynamics simulation of hen egg white lysozyme and normal-mode analyses, we find vibrations with large overlap with functional displacements are surprisingly concentrated in narrow frequency bands. These bands are not apparent in either the ensemble averaged vibrational density of states or isotropic absorption. However, in the case of the ensemble averaged anisotropic absorption, there is persistent spectral structure and overlap between this structure and the functional displacement frequency bands. We systematically lay out heuristics for calculating the spectra robustly, including the need for statistical sampling of the protein and inclusion of adequate water in the spectral calculation. The results show the congested spectrum of these complex molecules obscures important frequency bands associated with function and reveal a method to overcome this congestion by combining structurally sensitive spectroscopy with robust normal mode ensemble analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

CHARMM-GUI Free Energy Calculator for Absolute and Relative Ligand Solvation and Binding Free Energy Simulations

Alchemical free energy simulations have long been utilized to predict free energy changes for binding affinity and solubility of small molecules. However, while the theoretical foundation of these methods is well established, seamlessly handling many of the practical aspects regarding the preparation of the different thermodynamic end states of complex molecular systems and the numerous processing scripts often remains a burden for successful applications. Here, we present CHARMM-GUI Free Energy Calculator (http://www.charmm-gui.org/input/fec) that provides various alchemical free energy perturbation molecular dynamics (FEP/MD) systems with input and post-processing scripts for NAMD and GENESIS. Four submodules are available: Absolute Ligand Binder (for absolute ligand binding FEP/MD), Relative Ligand Binder (for relative ligand binding FEP/MD), Absolute Ligand Solvator (for absolute ligand solvation FEP/MD), and Relative Ligand Solvator (for relative ligand solvation FEP/MD). Each module is designed to build multiple systems of a set of selected ligands at once for high-throughput FEP/MD simulations. The capability of Free Energy Calculator is illustrated by absolute and relative solvation FEP/MD of a set of ligands and absolute and relative binding FEP/MD of a set of ligands for T4-lysozyme in solution and the adenosine A 2A receptor in a membrane. The calculated free energy values are overall consistent with the experimental and published free energy results (within ~1 kcal/mol). We hope that Free Energy Calculator is useful to carry out high-throughput FEP/MD simulations in the field of biomolecular sciences and drug discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Nuclear–Electronic Orbital Quantum Mechanical/Molecular Mechanical Real-Time Dynamics

Simulating the nuclear–electronic quantum dynamics of large-scale molecular systems in the condensed phase is key for studying biologically and chemically important processes such as proton transfer and proton-coupled electron transfer reactions. Herein, the real-time nuclear–electronic orbital time-dependent density functional theory (RT-NEO-TDDFT) approach is combined with a hybrid quantum mechanical/molecular mechanical (QM/MM) strategy to enable the accurate description of coupled nuclear–electronic quantum dynamics in the presence of heterogeneous environments such as solvent or proteins. The densities of the electrons and quantum protons are propagated in real time, while the other nuclei are propagated classically on the instantaneous electron–proton vibronic surface. This approach is applied to phenol bound to lysozyme, intramolecular proton transfer in malonaldehyde, and nonequilibrium excited-state intramolecular proton transfer in o-hydroxybenzaldehyde. Furthermore, these examples illustrate that the RT-NEO-TDDFT framework, coupled with an atomistic representation of the environment, allows the simulation of condensed-phase systems that exhibit significant nuclear quantum effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗