Investigating substitutions in antibody–antigen complexes using molecular dynamics: a case study with broad-spectrum, Influenza A antibodies

In studying the binding of host antibodies to the surface antigens of pathogens, the structural and functional characterization of antibody–antigen complexes by X-ray crystallography and binding assay is important. However, the characterization requires experiments that are typically time consuming and expensive: thus, many antibody–antigen complexes are under-characterized. For vaccine development and disease surveillance, it is often vital to assess the impact of amino acid substitutions on antibody binding. For example, are there antibody substitutions capable of improving binding without a loss of breadth, or antigen substitutions that lead to antigenic escape? The questions cannot be answered reliably from sequence variation alone, exhaustive substitution assays are usually impractical, and alanine scans provide at best an incomplete identification of the critical residue–residue interactions. Here, we show that, given an initial structure of an antibody bound to an antigen, molecular dynamics simulations using the energy method molecular mechanics with Generalized Born surface area (MM/GBSA) can model the impact of single amino acid substitutions on antibody–antigen binding energy. We apply the technique to three broad-spectrum antibodies to influenza A hemagglutinin and examine both previously characterized and novel variant strains observed in the human population that may give rise to antigenic escape. We find that in some cases the impact of a substitution is local, while in others it causes a reorientation of the antibody with wide-ranging impact on residue–residue interactions: this explains, in part, why the change in chemical properties of a residue can be, on its own, a poor predictor of overall change in binding energy. Our estimates are in good agreement with experimental results—indeed, they approximate the degree of agreement between different experimental techniques. Simulations were performed on commodity computer hardware; hence, this approach has the potential to be widely adopted by those undertaking infectious disease research. Novel aspects of this research include the application of MM/GBSA to investigate binding between broadly binding antibodies and a viral glycoprotein; the development of an approach for visualizing substrate–ligand interactions; and the use of experimental assay data to rescale our predictions, allowing us to make inferences about absolute, as well as relative, changes in binding energy

Ozone database in support of CMIP5 simulations: results and corresponding radiative forcing

A continuous tropospheric and stratospheric vertically resolved ozone time series, from 1850 to 2099, has been generated to be used as forcing in global climate models that do not include interactive chemistry. A multiple linear regression analysis of SAGE I+II satellite observations and polar ozonesonde measurements is used for the stratospheric zonal mean dataset during the well-observed period from 1979 to 2009. In addition to terms describing the mean annual cycle, the regression includes terms representing equivalent effective stratospheric chlorine (EESC) and the 11-yr solar cycle variability. The EESC regression fit coefficients, together with pre-1979 EESC values, are used to extrapolate the stratospheric ozone time series backward to 1850. While a similar procedure could be used to extrapolate into the future, coupled chemistry climate model (CCM) simulations indicate that future stratospheric ozone abundances are likely to be significantly affected by climate change, and capturing such effects through a regression model approach is not feasible. Therefore, the stratospheric ozone dataset is extended into the future (merged in 2009) with multimodel mean projections from 13 CCMs that performed a simulation until 2099 under the SRES (Special Report on Emission Scenarios) A1B greenhouse gas scenario and the A1 adjusted halogen scenario in the second round of the Chemistry-Climate Model Validation (CCMVal-2) Activity. The stratospheric zonal mean ozone time series is merged with a three-dimensional tropospheric data set extracted from simulations of the past by two CCMs (CAM3.5 and GISSPUCCINI)and of the future by one CCM (CAM3.5). The future tropospheric ozone time series continues the historical CAM3.5 simulation until 2099 following the four different Representative Concentration Pathways (RCPs). Generally good agreement is found between the historical segment of the ozone database and satellite observations, although it should be noted that total column ozone is overestimated in the southern polar latitudes during spring and tropospheric column ozone is slightly underestimated. Vertical profiles of tropospheric ozone are broadly consistent with ozonesondes and in-situ measurements, with some deviations in regions of biomass burning. The tropospheric ozone radiative forcing (RF) from the 1850s to the 2000s is 0.23Wm−2, lower than previous results. The lower value is mainly due to (i) a smaller increase in biomass burning emissions; (ii) a larger influence of stratospheric ozone depletion on upper tropospheric ozone at high southern latitudes; and possibly (iii) a larger influence of clouds (which act to reduce the net forcing) compared to previous radiative forcing calculations. Over the same period, decreases in stratospheric ozone, mainly at high latitudes, produce a RF of −0.08Wm−2, which is more negative than the central Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4) value of −0.05Wm−2, but which is within the stated range of −0.15 to +0.05Wm−2. The more negative value is explained by the fact that the regression model simulates significant ozone depletion prior to 1979, in line with the increase in EESC and as confirmed by CCMs, while the AR4 assumed no change in stratospheric RF prior to 1979. A negative RF of similar magnitude persists into the future, although its location shifts from high latitudes to the tropics. This shift is due to increases in polar stratospheric ozone, but decreases in tropical lower stratospheric ozone, related to a strengthening of the Brewer-Dobson circulation, particularly through the latter half of the 21st century. Differences in trends in tropospheric ozone among the four RCPs are mainly driven by different methane concentrations, resulting in a range of tropospheric ozone RFs between 0.4 and 0.1Wm−2 by 2100. The ozone dataset described here has been released for the Coupled Model Intercomparison Project (CMIP5) model simulations in netCDF Climate and Forecast (CF) Metadata Convention at the PCMDI website (http://cmip-pcmdi.llnl.gov/)

Laser operation of an Nd:Gd₃Ga₅O₁₂ thin-film optical waveguide fabricated by pulsed laser deposition

We report the laser operation of a thin-film waveguide structure grown by the pulsed laser deposition technique. A 2.7-µm-thick crystalline film of neodymium doped Gd3Ga5O12 (Nd:GGG) lases at a wavelength centered at 1.06µm when pumped by a Ti:sapphire laser at 808 nm

The power of quantum systems on a line

We study the computational strength of quantum particles (each of finite dimensionality) arranged on a line. First, we prove that it is possible to perform universal adiabatic quantum computation using a one-dimensional quantum system (with 9 states per particle). This might have practical implications for experimentalists interested in constructing an adiabatic quantum computer. Building on the same construction, but with some additional technical effort and 12 states per particle, we show that the problem of approximating the ground state energy of a system composed of a line of quantum particles is QMA-complete; QMA is a quantum analogue of NP. This is in striking contrast to the fact that the analogous classical problem, namely, one-dimensional MAX-2-SAT with nearest neighbor constraints, is in P. The proof of the QMA-completeness result requires an additional idea beyond the usual techniques in the area: Not all illegal configurations can be ruled out by local checks, so instead we rule out such illegal configurations because they would, in the future, evolve into a state which can be seen locally to be illegal. Our construction implies (assuming the quantum Church-Turing thesis and that quantum computers cannot efficiently solve QMA-complete problems) that there are one-dimensional systems which take an exponential time to relax to their ground states at any temperature, making them candidates for being one-dimensional spin glasses.Comment: 21 pages. v2 has numerous corrections and clarifications, and most importantly a new author, merged from arXiv:0705.4067. v3 is the published version, with additional clarifications, publisher's version available at http://www.springerlink.co

Star and Planet Formation with ALMA: an Overview

Submillimeter observations with ALMA will be the essential next step in our understanding of how stars and planets form. Key projects range from detailed imaging of the collapse of pre-stellar cores and measuring the accretion rate of matter onto deeply embedded protostars, to unravelling the chemistry and dynamics of high-mass star-forming clusters and high-spatial resolution studies of protoplanetary disks down to the 1 AU scale.Comment: Invited review, 8 pages, 5 figures; to appear in the proceedings of "Science with ALMA: a New Era for Astrophysics". Astrophysics & Space Science, in pres

(Sub)mm Interferometry Applications in Star Formation Research

This contribution gives an overview about various applications of (sub)mm interferometry in star formation research. The topics covered are molecular outflows, accretion disks, fragmentation and chemical properties of low- and high-mass star-forming regions. A short outlook on the capabilities of ALMA is given as well.Comment: 20 pages, 7 figures, in proceedings to "2nd European School on Jets from Young Star: High Angular Resolution Observations". A high-resolution version of the paper can be found at http://www.mpia.de/homes/beuther/papers.htm

Extending Torelli map to toroidal compactifications of Siegel space

It has been known since the 1970s that the Torelli map $M_g \to A_g$, associating to a smooth curve its jacobian, extends to a regular map from the Deligne-Mumford compactification $\bar{M}_g$ to the 2nd Voronoi compactification $\bar{A}_g^{vor}$. We prove that the extended Torelli map to the perfect cone (1st Voronoi) compactification $\bar{A}_g^{perf}$ is also regular, and moreover $\bar{A}_g^{vor}$ and $\bar{A}_g^{perf}$ share a common Zariski open neighborhood of the image of $\bar{M}_g$. We also show that the map to the Igusa monoidal transform (central cone compactification) is NOT regular for $g\ge9$; this disproves a 1973 conjecture of Namikawa.Comment: To appear in Inventiones Mathematica

First Measurement of Z/gamma* Production in Compton Scattering of Quasi-real Photons

We report the first observation of Z/gamma* production in Compton scattering of quasi-real photons. This is a subprocess of the reaction e+e- to e+e-Z/gamma*, where one of the final state electrons is undetected. Approximately 55 pb-1 of data collected in the year 1997 at an e+e- centre-of-mass energy of 183 GeV with the OPAL detector at LEP have been analysed. The Z/gamma* from Compton scattering has been detected in the hadronic decay channel. Within well defined kinematic bounds, we measure the product of cross-section and Z/gamma* branching ratio to hadrons to be (0.9+-0.3+-0.1) pb for events with a hadronic mass larger than 60 GeV, dominated by (e)eZ production. In the hadronic mass region between 5 GeV and 60 GeV, dominated by (e)egamma* production, this product is found to be (4.1+-1.6+-0.6) pb. Our results agree with the predictions of two Monte Carlo event generators, grc4f and PYTHIA.Comment: 18 pages, LaTeX, 5 eps figures included, submitted to Physics Letters

Search for Higgs Bosons in e+e- Collisions at 183 GeV

The data collected by the OPAL experiment at sqrts=183 GeV were used to search for Higgs bosons which are predicted by the Standard Model and various extensions, such as general models with two Higgs field doublets and the Minimal Supersymmetric Standard Model (MSSM). The data correspond to an integrated luminosity of approximately 54pb-1. None of the searches for neutral and charged Higgs bosons have revealed an excess of events beyond the expected background. This negative outcome, in combination with similar results from searches at lower energies, leads to new limits for the Higgs boson masses and other model parameters. In particular, the 95% confidence level lower limit for the mass of the Standard Model Higgs boson is 88.3 GeV. Charged Higgs bosons can be excluded for masses up to 59.5 GeV. In the MSSM, mh > 70.5 GeV and mA > 72.0 GeV are obtained for tan{beta}>1, no and maximal scalar top mixing and soft SUSY-breaking masses of 1 TeV. The range 0.8 < tanb < 1.9 is excluded for minimal scalar top mixing and m{top} < 175 GeV. More general scans of the MSSM parameter space are also considered.Comment: 49 pages. LaTeX, including 33 eps figures, submitted to European Physical Journal

A Measurement of the Product Branching Ratio f(b->Lambda_b).BR(Lambda_b->Lambda X) in Z0 Decays

The product branching ratio, f(b->Lambda_b).BR(Lambda_b->Lambda X), where Lambda_b denotes any weakly-decaying b-baryon, has been measured using the OPAL detector at LEP. Lambda_b are selected by the presence of energetic Lambda particles in bottom events tagged by the presence of displaced secondary vertices. A fit to the momenta of the Lambda particles separates signal from B meson and fragmentation backgrounds. The measured product branching ratio is f(b->Lambda_b).BR(Lambda_b->Lambda X) = (2.67+-0.38(stat)+0.67-0.60(sys))% Combined with a previous OPAL measurement, one obtains f(b->Lambda_b).BR(Lambda_b->Lambda X) = (3.50+-0.32(stat)+-0.35(sys))%.Comment: 16 pages, LaTeX, 3 eps figs included, submitted to the European Physical Journal