A Case of an Interventricular Septum Pseudoaneurysm With Perforation Mimicking a Ventricular Septal Defect.

Ventricular pseudoaneurysm (PSA) is a ventricular outpouching contained by adherent pericardium or myocardial scar tissue and represents a rare but potentially fatal complication of acute myocardial infarction (AMI). The vast majority of cases involve the left ventricular apex, in the area of infarct. It is extremely rare to see PSA formation within the interventricular septum (IVS).

Unusual Presentation of Spinal Myoclonus After Neuraxial Anesthesia in a Patient With Multiple Sclerosis: A Case Report.

This case report highlights a unique instance of spinal myoclonus after neuraxial anesthesia. It aims to inform anesthesiology providers, enhancing their ability to identify, manage, and potentially prevent similar outcomes in patients at risk.

Molecular modeling of nucleic acid structure: electrostatics and solvation.

This unit presents an overview of computer simulation techniques as applied to nucleic acid systems, ranging from simple in vacuo molecular modeling techniques to more complete all-atom molecular dynamics treatments that include an explicit representation of the environment. The third in a series of four units, this unit focuses on critical issues in solvation and the treatment of electrostatics.

Molecular modeling of nucleic acid structure: energy and sampling.

An overview of computer simulation techniques as applied to nucleic acid systems is presented. This unit expands an accompanying overview unit (UNIT 7.5) by discussing methods used to treat the energy and sample representative configurations.

Molecular modeling of nucleic acid structure.

This unit is the first in a series of four units covering the analysis of nucleic acid structure by molecular modeling. This unit provides an overview of computer simulation of nucleic acids. Topics include the static structure model, computational graphics and energy models, generation of an initial model, and characterization of the overall three-dimensional structure.

Molecular modeling of nucleic acid structure: setup and analysis.

The last in a set of units by these authors, this unit addresses some important remaining questions about molecular modeling of nucleic acids. It describes how to choose an appropriate molecular mechanics force field; how to set up and equilibrate the system for accurate simulation of a nucleic acid in an explicit solvent by molecular dynamics or Monte Carlo simulation; and how to analyze molecular dynamics trajectories.

New-generation amber united-atom force field.

We have developed a new-generation Amber united-atom force field for simulations involving highly demanding conformational sampling such as protein folding and protein-protein binding. In the new united-atom force field, all hydrogens on aliphatic carbons in all amino acids are united with carbons except those on Calpha. Our choice of explicit representation of all protein backbone atoms aims at minimizing perturbation to protein backbone conformational distributions and to simplify development of backbone torsion terms.

Chelate effect in cyclodextrin dimers: a computational (MD, MM/PBSA, and MM/GBSA) study.

The complexation of an adamantyl-phosphate derivative with one beta-cyclodextrin, with two beta-cyclodextrins, and with two beta-cyclodextrins dimerized with a disulfide bridge was studied by computational methods (MD, MM/PBSA, and MM/GBSA) to analyze and rationalize the chelate effect. Although this effect is usually explained by invoking favorable entropy contribution due to the preorganization of the ligand, it has been determined experimentally that in this case it is enthalpy-driven. The computational results are in accord with this finding, although the entropy contribution due to the solvent structural organization around the complex is crucial for the final estimates of the free energy of complexation.

Automatic atom type and bond type perception in molecular mechanical calculations.

In molecular mechanics (MM) studies, atom types and/or bond types of molecules are needed to determine prior to energy calculations. We present here an automatic algorithm of perceiving atom types that are defined in a description table, and an automatic algorithm of assigning bond types just based on atomic connectivity. The algorithms have been implemented in a new module of the AMBER packages.

The open nucleotide pocket of the profilin/actin x-ray structure is unstable and closes in the absence of profilin.

The open nucleotide pocket conformation of actin in the profilin:actinCaATP x-ray structure has been hypothesized to be a crucial intermediate for nucleotide exchange in the actin depolymerization/polymerization cycle. The requirement for ancillary modification of actin for crystallization leads to ambiguities in this interpretation, however. We have used molecular dynamics simulations to model the thermodynamic properties of the actin x-ray structure, outside the crystal lattice, in an aqueous environment with profilin removed.

Hierarchical database screenings for HIV-1 reverse transcriptase using a pharmacophore model, rigid docking, solvation docking, and MM-PB/SA.

In this work, an efficient strategy was presented to search drug leads for human immunodeficiency virus type 1 reverse transcriptase (HIV-1 RT) using hierarchical database screenings, which included a pharmacophore model, multiple-conformation rigid docking, solvation docking, and molecular mechanics-Poisson-Boltzmann/surface area (MM-PB/SA) sequentially. Encouraging results were achieved in searching a refined available chemical directory (ACD) database: the enrichment factor after the first three filters was estimated to be 25-fold; the hit rate for all the four filters was predicted to be 41% in a control test using 37 known HIV-1 non-nucleoside reverse transcriptase inhibitors; 10 out of 30 promising solvation-docking hits had MM-PB/SA binding free energies better than -6.8 kcal/mol and the best one, HIT15, had -17.

Prediction of pKa shifts in proteins using a combination of molecular mechanical and continuum solvent calculations.

The prediction of pKa shifts of ionizable groups in proteins is of great relevance for a number of important biological phenomena. We present an implementation of the MM-GBSA approach, which combines molecular mechanical (MM) and generalized Born (GB) continuum solvent energy terms, to the calculation of pKa values of a panel of nine proteins, including 69 individual comparisons with experiment. While applied so far mainly to the calculation of biomolecular binding free energies, we show that this method can also be used for the estimation of protein pKa shifts, with an accuracy around 1 pKa unit, even for strongly shifted residues.

Development and testing of a general amber force field.

We describe here a general Amber force field (GAFF) for organic molecules. GAFF is designed to be compatible with existing Amber force fields for proteins and nucleic acids, and has parameters for most organic and pharmaceutical molecules that are composed of H, C, N, O, S, P, and halogens. It uses a simple functional form and a limited number of atom types, but incorporates both empirical and heuristic models to estimate force constants and partial atomic charges.

Enhanced ab initio protein folding simulations in Poisson-Boltzmann molecular dynamics with self-guiding forces.

We have investigated the sampling efficiency in molecular dynamics with the PB implicit solvent when self-guiding forces are added. Compared with a high-temperature dynamics simulation, the use of self-guiding forces in room-temperature dynamics is found to be rather efficient as measured by potential energy fluctuation, gyration radius fluctuation, backbone RMSD fluctuation, number of unique clusters, and distribution of low RMSD structures over simulation time. Based on the enhanced sampling method, we have performed ab initio folding simulations of two small proteins, betabetaalpha1 and villin headpiece.

Investigation of neuraminidase-substrate recognition using molecular dynamics and free energy calculations.

Development of the new generation of therapeutics against the influenza viral coat protein neuraminidase is a response to the continuing threat of influenza epidemics. A variety of structurally similar compounds have been reported that vary greatly in their ability to inhibit neuraminidase, a critical enzyme that cleaves sialic acid and promotes virion release. To determine how neuraminidase exhibits this wide range of affinities with structurally similar compounds, molecular dynamic simulations, coupled with free energy calculations, were used to determine the binding components of a series of neuraminidase inhibitors.

A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations.

Molecular mechanics models have been applied extensively to study the dynamics of proteins and nucleic acids. Here we report the development of a third-generation point-charge all-atom force field for proteins. Following the earlier approach of Cornell et al.

Direct hydroxide attack is a plausible mechanism for amidase antibody 43C9.

Direct hydroxide attack on the scissile carbonyl of the substrate has been suggested as a likely mechanism for esterase antibodies elicited by phosphonate haptens, which mimic the transition states for the alkaline hydrolysis of esters.1 The unique amidase activity of esterase antibody 43C9 has been attributed to nucleophilic attack by an active-site histidine residue.2 Yet, the active site of 43C9 is strikingly similar to those of other esterase antibodies, particularly 17E8.

Closing of the nucleotide pocket of kinesin-family motors upon binding to microtubules.

We have used adenosine diphosphate analogs containing electron paramagnetic resonance (EPR) spin moieties and EPR spectroscopy to show that the nucleotide-binding site of kinesin-family motors closes when the motor.diphosphate complex binds to microtubules. Structural analyses demonstrate that a domain movement in the switch 1 region at the nucleotide site, homologous to domain movements in the switch 1 region in the G proteins [heterotrimeric guanine nucleotide-binding proteins], explains the EPR data.

Free energy calculations for theophylline binding to an RNA aptamer: Comparison of MM-PBSA and thermodynamic integration methods.

We have applied the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method (J. Srinivasan, T. E.

The intramolecular mechanism for the second proton transfer in triosephosphate isomerase (TIM): a QM/FE approach.

The intramolecular mechanism we earlier proposed [Alagona, G.; Desmeules, P.; Ghio, C.

pKa, MM, and QM studies of mechanisms of beta-lactamases and penicillin-binding proteins: acylation step.

The acylation step of the catalytic mechanism of beta-lactamases and penicillin-binding proteins (PBPs) has been studied with various approaches. The methods applied range from molecular dynamics (MD) simulations to multiple titration calculations using the Poisson-Boltzmann approach to quantum mechanical (QM) methods. The mechanism of class A beta-lactamases was investigated in the greatest detail.

Molecular dynamics simulation study of the negative correlation in antibody AZ28-catalyzed oxy-cope rearrangement.

The oxy-Cope rearrangement reaction in the antibody AZ28 is investigated using ab initio molecular orbital calculations and molecular mechanical molecular dynamics simulations. This antibody, AZ28, is known as one of the few systems where the mature catalytic antibody shows a negative correlation between the transition state analogue (TSA) binding affinity and the catalytic rate of the oxy-Cope rearrangement compared to the germ line catalytic antibody. The ab initio optimized structure shows that the transition state structure has a more planar configuration than the TSA.

Simulating proteins at constant pH: An approach combining molecular dynamics and Monte Carlo simulation.

For the structure and function of proteins, the pH of the solution is one of the determining parameters. Current molecular dynamics (MD) simulations account for the solution pH only in a limited way by keeping each titratable site in a chosen protonation state. We present an algorithm that generates trajectories at a Boltzmann distributed ensemble of protonation states by a combination of MD and Monte Carlo (MC) simulation.

Conformational preferences of substituted prolines in the collagen triple helix.

Researchers have recently questioned the role hydroxylated prolines play in stabilizing the collagen triple helix. To address these issues, we have developed new molecular mechanics parameters for the simulation of peptides containing 4(R)-fluoroproline (Flp), 4(R)-hydroxyproline (Hyp), and 4(R)-aminoproline (Amp). Simulations of peptides based on these parameters can be used to determine the components that stabilize hydroxyproline over proline in the triple helix.