Neutralizing antibodies (Nab) are a principal component of an effective human immune response to many pathogens, yet their role in HIV-1 infection is unclear. To gain a better understanding of this role, we examined plasma from patients with acute HIV infection. Here we report the detection of autologous Nab as early as 52 days after detection of HIV-specific antibodies. The viral inhibitory activity of Nab resulted in complete replacement of neutralization-sensitive virus by successive populations of resistant virus. Escape virus contained mutations in the env gene that were unexpectedly sparse, did not map generally to known neutralization epitopes, and involved primarily changes in N-linked glycosylation. This pattern of escape, and the exceptional density of HIV-1 envelope glycosylation generally, led us to postulate an evolving 'glycan shield' mechanism of neutralization escape whereby selected changes in glycan packing prevent Nab binding but not receptor binding. Direct support for this model was obtained by mutational substitution showing that Nab-selected alterations in glycosylation conferred escape from both autologous antibody and epitope-specific monoclonal antibodies. The evolving glycan shield thus represents a new mechanism contributing to HIV-1 persistence in the face of an evolving antibody repertoire.

Antibody neutralization and escape by HIV-1

We present an implementation of generalized Born implicit solvent all-atom classical molecular dynamics (MD) within the AMBER program package that runs entirely on CUDA enabled NVIDIA graphics processing units (GPUs). We discuss the algorithms that are used to exploit the processing power of the GPUs and show the performance that can be achieved in comparison to simulations on conventional CPU clusters. The implementation supports three different precision models in which the contributions to the forces are calculated in single precision floating point arithmetic but accumulated in double precision (SPDP), or everything is computed in single precision (SPSP) or double precision (DPDP). In addition to performance, we have focused on understanding the implications of the different precision models on the outcome of implicit solvent MD simulations. We show results for a range of tests including the accuracy of single point force evaluations and energy conservation as well as structural properties pertainining to protein dynamics. The numerical noise due to rounding errors within the SPSP precision model is sufficiently large to lead to an accumulation of errors which can result in unphysical trajectories for long time scale simulations. We recommend the use of the mixed-precision SPDP model since the numerical results obtained are comparable with those of the full double precision DPDP model and the reference double precision CPU implementation but at significantly reduced computational cost. Our implementation provides performance for GB simulations on a single desktop that is on par with, and in some cases exceeds, that of traditional supercomputers.

/pdf/routine-microsecond-molecular-dynamics-simulations-with-41c2fb6wpv.pdf

Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 1. Generalized Born

We study the system-size dependence of translational diffusion coefficients and viscosities in molecular dynamics simulations under periodic boundary conditions. Simulations of water under ambient conditions and a Lennard-Jones (LJ) fluid show that the diffusion coefficients increase strongly as the system size increases. We test a simple analytic correction for the system-size effects that is based on hydrodynamic arguments. This correction scales as N-1/3, where N is the number of particles. For a cubic simulation box of length L, the diffusion coefficient corrected for system-size effects is D0 = DPBC + 2.837297kBT/(6πηL), where DPBC is the diffusion coefficient calculated in the simulation, kB the Boltzmann constant, T the absolute temperature, and η the shear viscosity of the solvent. For water, LJ fluids, and hard-sphere fluids, this correction quantitatively accounts for the system-size dependence of the calculated self-diffusion coefficients. In contrast to diffusion coefficients, the shear viscos...

System-Size Dependence of Diffusion Coefficients and Viscosities from Molecular Dynamics Simulations with Periodic Boundary Conditions

The AMBER lipid force field has been updated to create Lipid14, allowing tensionless simulation of a number of lipid types with the AMBER MD package. The modular nature of this force field allows numerous combinations of head and tail groups to create different lipid types, enabling the easy insertion of new lipid species. The Lennard-Jones and torsion parameters of both the head and tail groups have been revised and updated partial charges calculated. The force field has been validated by simulating bilayers of six different lipid types for a total of 0.5 μs each without applying a surface tension; with favorable comparison to experiment for properties such as area per lipid, volume per lipid, bilayer thickness, NMR order parameters, scattering data, and lipid lateral diffusion. As the derivation of this force field is consistent with the AMBER development philosophy, Lipid14 is compatible with the AMBER protein, nucleic acid, carbohydrate, and small molecule force fields.

/pdf/lipid14-the-amber-lipid-force-field-1gj0t3c4bx.pdf

Lipid14: The Amber Lipid Force Field

Obtaining the correct balance of secondary structure propensities is a central priority in protein force-field development. Given that current force fields differ significantly in their alpha-helical propensities, a correction to match experimental results would be highly desirable. We have determined simple backbone energy corrections for two force fields to reproduce the fraction of helix measured in short peptides at 300 K. As validation, we show that the optimized force fields produce results in excellent agreement with nuclear magnetic resonance experiments for folded proteins and short peptides not used in the optimization. However, despite the agreement at ambient conditions, the dependence of the helix content on temperature is too weak, a problem shared with other force fields. A fit of the Lifson-Roig helix-coil theory shows that both the enthalpy and entropy of helix formation are too small: the helix extension parameter w agrees well with experiment, but its entropic and enthalpic components are both only about half the respective experimental estimates. Our structural and thermodynamic analyses point toward the physical origins of these shortcomings in current force fields, and suggest ways to address them in future force-field development.

/pdf/optimized-molecular-dynamics-force-fields-applied-to-the-1g13okgc3h.pdf

Optimized molecular dynamics force fields applied to the helix-coil transition of polypeptides.

End-to-end contact formation rates of several peptides were recently measured by tryptophan triplet quenching (Lapidus et al. Proc. Natl. Acad. Sci. U.S.A. 2000, 97, 7220). Motivated by these experiments, we study loop-closure kinetics for two peptides of different lengths, Cys-(Ala-Gly-Gln)n-Trp (n = 1, 2), in multiple all-atom explicit-solvent molecular dynamics simulations with different initial conditions and force fields. In 150 simulations of approximately 20 ns each, we collect data covering 1.0 and 0.8 μs for the penta-peptide simulated with the AMBER and CHARMM force fields, respectively, and about 0.5 μs each with the two force fields for the octa-peptide. These extensive simulations allow us to analyze the dynamics of peptides in the unfolded state with atomic resolution, thus probing early events in protein folding, and to compare molecular dynamics simulations directly with experiment. The calculated lifetimes of the tryptophan triplet state are in the range of 50−100 ns, in agreement with ex...

https://faculty.uml.edu/vbarsegov/teaching/bioinformatics/papers/yeh.pdf

Peptide loop-closure kinetics from microsecond molecular dynamics simulations in explicit solvent.

Diffusion and electrophoretic mobility of single-stranded RNA from molecular dynamics simulations.

We analyzed the thermodynamic and structural determinants of indolicidin interactions with eukaryotic and prokaryotic cell membranes using a series of atomistically detailed molecular dynamics simulations. We used quartz-supported bilayers with two different compositions of zwitterionic and anionic phospholipids as model eukaryotic and prokaryotic cell membranes. Indolicidin was preferentially attracted to the model prokaryotic cell membrane in contrast to the weak adsorption on the eukaryotic membrane. The nature of the indolicidin surface adsorption depended on an electrostatic guiding component, an attractive enthalpic component derived from van der Waals interactions, and a balance between entropic factors related to peptide confinement at the interface and counterion release from the bilayer surface. Thus, whereas we attributed the specificity of the indolicidin/membrane interaction to electrostatics, these interactions were not the sole contributors to the free energy of adsorption. Instead, a balance between an attractive van der Waals enthalpic component and a repulsive entropic component determined the overall strength of indolicidin adsorption.

In-Chul Yeh

Papers

System-Size Dependence of Diffusion Coefficients and Viscosities from Molecular Dynamics Simulations with Periodic Boundary Conditions

Peptide loop-closure kinetics from microsecond molecular dynamics simulations in explicit solvent.

Diffusion and electrophoretic mobility of single-stranded RNA from molecular dynamics simulations.

Free Energy Difference in Indolicidin Attraction to Eukaryotic and Prokaryotic Model Cell Membranes B

Molecular Dynamics Simulation of Solid-Supported Lipid Bilayers