PubMed Health⌕ Search

Biomedical subjects

Edward N Brothers

Publications and source records attributed to Edward N Brothers.

8 recordsLinked to original sources

Linear-scaling calculation of static and dynamic polarizabilities in Hartree-Fock and density functional theory for periodic systems.

We present a linear-scaling method for analytically calculating static and dynamic polarizabilities with Hartree-Fock and density functional theory, using Gaussian orbitals and periodic boundary conditions. Our approach uses the direct space fast multipole method to evaluate the long-range Coulomb contributions. For exact exchange, we use efficient screening techniques developed for energy calculations. We then demonstrate the capabilities of our approach with benchmark calculations on one-, two-, and three-dimensional systems.

Journal Article↗

Longitudinal polarizability of carbon nanotubes.

The longitudinal polarizabilities of carbon nanotubes are determined using first principles density functional theory. These results demonstrate that the polarizability per atom of a nanotube in the axial direction is primarily determined by the band gap. In fact, polarizability per atom versus inverse band gap yields a linear trend for all nanotubes and methods utilized in this study, creating a universal relationship for longitudinal polarizability. This can be explained by examining the terms in the sum over states equation used to determine polarizability and noting that the vast majority of the polarizability arises from a few elements near the band gap. This universal trend is then used with experimentally determined band gaps to predict the experimental polarizability of carbon nanotubes.

Journal Article↗

Coaxial carbon nanotubes as shielded nanowires.

We report transverse polarizabilities of coaxial carbon nanotubes using first principles density functional theory. These results demonstrate a shielding of the inner nanotube from electric fields by the outer nanotube. This study has implications for nanoelectronics, specifically for the possibility of using coaxial nanotubes as nanoelectrical wires. Shielding is predicted to be on the order of 95% by high-level polarizability calculations. This shielding occurs regardless of whether the outer nanotube is metallic or semiconducting. In addition, a series of calculations on coaxial nanotubes where the inner nanotube is not centered show that the shielding still occurs with approximately the same magnitude. These calculations therefore indicate that it would be possible to use a coaxial carbon nanotube as a shielded nanowire.

Journal Article↗

Hybrid QM/MM and DFT investigations of the catalytic mechanism and inhibition of the dinuclear zinc metallo-beta-lactamase CcrA from Bacteroides fragilis.

Based on hybrid QM/MM molecular dynamics simulation and density functional theoretical (DFT) calculations, we investigate the mechanistic and energetic features of the catalytic action of dizinc metallo-beta-lactamase CcrA from Bacteroides fragilis. The 200 ps QM/MM simulation of the CcrA enzyme in complex with nitrocefin shows that the substrate beta-lactam moiety is directed toward the active site dizinc center through the interactions of aminocarbonyl and carboxylate groups with the two active site zinc ions and the two conserved residues, Lys167 and Asn176. From the determination of the potential energy profile of a relevant enzymatic reaction model, it is found that the nucleophilic displacement reaction step proceeds with a low-barrier height, leading to the formation of an energetically favored reaction intermediate. The results also show that the high catalytic activity of the CcrA enzyme stems from a simultaneous operation of three catalytic components: activation of the bridging hydroxide nucleophile by zinc-coordinated Asp86; polarization of the substrate aminocarbonyl group by the first zinc ion; stabilization of the negative charge developed on the departing amide nitrogen by the second zinc ion. Consistent with the previous experimental finding that the proton-transfer reaction step is rate-limiting, the activation energy of the second step is found to be 1.6 kcal/mol higher than that of the first step. Finally, through an examination of the structural and energetic features of binding of a thiazolidinecarboxylic acid inhibitor to the active site dizinc center, a two-step inhibition mechanism involving a protonation-induced ligand exchange reaction is proposed for the inhibitory action of a tight-binding inhibitor possessing a thiol group.

Bacterial Proteins↗

PM3-compatible zinc parameters optimized for metalloenzyme active sites.

Recent studies have shown that semiempirical methods (e.g., PM3 and AM1) for zinc-containing compounds are unreliable for modeling structures containing zinc ions with ligand environments similar to those observed in zinc metalloenzymes. To correct these deficiencies a reparameterization of zinc at the PM3 level was undertaken. In this effort we included frequency corrected B3LYP/6-311G* zinc metalloenzyme ligand environments along with previously utilized experimental data. Average errors for the heats of formation have been reduced from 46.9 kcal/mol (PM3) to 14.2 kcal/mol for this new parameter set, termed ZnB for "Zinc, Biological." In addition, the new parameter sets predict geometries for the Bacillus fragilis active site model and other zinc metalloenzyme mimics that are qualitatively in agreement with high-level ab initio results, something existing parameter sets failed to do.

Bacillus↗

A finite temperature linear tetrahedron method for electronic structure calculations of periodic systems.

A finite-temperature linear tetrahedron method for electronic structure calculations of periodic systems is developed. When compared to widely used simple temperature broadening, the number of k points necessary for accurate integration at finite temperatures is reduced. The utility of the method is demonstrated with benchmark calculations on 1D, 2D, and 3D systems.

Journal Article↗

Fast semiempirical calculations for nuclear magnetic resonance chemical shifts: a divide-and-conquer approach.

A new approach to calculate nuclear magnetic resonance chemical shifts has been implemented at the semiempirical modified neglect of diatomic overlap level using gauge-including atomic orbitals. The perturbed density matrix with respect to the magnetic field is obtained by the diagonalization of the complex Fock matrix using the divide and conquer (DC) method, instead of by solving the computationally expensive coupled perturbed Hartree-Fock equations. Adopting the Patchkovskii and Thiel parameters [S. Patchkovskii and W. Thiel J. Comput. Chem. 20, 1220 (1999)], we were able to reproduce their results for small organic molecules. The errors introduced by DC method are negligible, as shown by the calculations on a series of polyalaine structures. Test calculations on proteins have demonstrated that our approach makes it possible to calculate chemical shifts routinely on systems with hundreds of atoms with good accuracy.

Journal Article↗

Insights into the structure and dynamics of the dinuclear zinc beta-lactamase site from Bacteroides fragilis.

Herein, we report quantum chemical calculations and molecular dynamics (MD) simulations of the dinuclear form of the Bacteroides fragilis zinc beta-lactamase. We studied four different configurations which differ in the protonation state of the Asp103 residue and in the presence or absence of a Zn1-OH-Zn2 bridge. The flexibility of the Zn1-OH-Zn2 bridge was studied by means of quantum mechanical (QM) calculations on cluster models while the relative stabilities of the different configurations were estimated from QM linear scaling calculations on the enzyme. Contacts between important residues (Cys104, Asp69, Lys185, etc.), the solvation of the zinc ions, and the conformation of the active site beta-hairpin loop were characterized by the MD analyses. The influence of the buried sodium ion close to the Zn2 position was investigated by carrying out a secondary simulation where the sodium ion was replaced with an internal water molecule. The comparative structural analyses among the different MD trajectories augmented with energetic calculations have demonstrated that the B. fragilis protein efficiently binds the internal Na(+) ion observed crystallographically. Moreover, we found that when Asp103 is unprotonated, a rigid Zn1-OH-Zn2 bridge results, while for neutral Asp103, a fluctuating Zn1-Zn2 distance was possible via the breaking and formation of the Zn1-OH-Zn2 bridge. The mechanistic implications of these observations are discussed in detail.

Aspartic Acid↗